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A short introduction is given on the functional renormalization group method, putting emphasis 
on its nonperturbative aspects. The method enables to find nontrivial fixed points in quantum field 
theoretic models which make them free from divergences and leads to the concept of asymptotic 
safety. It can be considered as a generalization of the asymptotic freedom which plays a key role 
in the perturbative renormalization. We summarize and give a short discussion of some important 
models, which are asymptotically safe such as the Gross-Neveu model, the nonlinear a model, the 
sine-Gordon model, and the model of quantum Einstein gravity. We also give a detailed analysis 
" of infrared behavior of the models where a spontaneous symmetry breaking takes place. The deep 

infrared behavior of the broken phase cannot be treated within the framework of perturbative calcu- 
lations. We demonstrate that there exists an infrared fixed point in the broken phase which creates 
£SJ ' a new scaling regime there, however its structure is hidden by the singularity of the renormalization 

group equations. The phase spaces of these models show several similar properties, namely the 
models has the same phase and fixed point structure. These results can only be uncovered by the 
functional renormalization group method. 
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I. INTRODUCTION 



The concept of asymptotic safety can be understood in the framework of quantum field theory, where we treat 
quantum systems with many degrees of freedom and many orders of magnitude in length, or in energy, momentum. 
Usually we have an assumption of describing the high energy/ultra violet (UV), or short distance/microscopic in- 
teractions among the elementary particles where a high degree of symmetry exists. However the measurements are 
typically performed at low energies, in the infrared (IR) regime, therefore we need the theory, where those quantum 
fluctuations or modes are taken into account which energy is between the UV and IR energy scales. The functional 
(/"■) . renormalization group (RG) method is one of the best candidate to take into account the quantum fluctuations step 
t-H ■ by step, one by one, systematically [1-7]. The method provides us a partial differential differential equations for the 
^sl" | action of the model. As an initial condition we need the UV action, and the solution of the equation gives us the IR 
■ one which can be identified with the effective potential. If these scales are close to each other, then even a perturbative 
treatment can be satisfactory. However, generally one cannot conclude to the IR behavior from the UV one. 

The effective action can be derived from the blocked action. We assume (especially in scalar models), that the 
blocked and the effective action has the same functional structure. We usually perform a field expansion to the 
effective action which results in many interacting terms, and each of them is multiplied by a coupling. The original 
partial differential equation can be converted to a system of ordinary differential equations w.r.t. the couplings. These 
equations are called evolution or flow equations. The initial conditions are the UV values of the couplings and the 
results of the equations are the IR ones. The RG method typically gives highly nonlinear evolution equations, giving 
flows which cannot be recovered by a perturbative RG treatment. 

Moreover the perturbative RG calculations are performed in the vicinity of the Gaussian fixed point (GFP). The 
GFP is the origin of the phase space which is spanned by the (dimensionless) couplings. The UV effective action is 
given by small values of the initial couplings, so the kinetic term dominates the UV physics. The UV scaling of the 
couplings are UV attractive or IR repulsive (they go away from the fixed point if we consider the direction of the RG 
'time' scale) if the model is perturbatively renormalizable. Otherwise it may happen that the GFP is a hyperbolic 
point which means that there are directions in the phase space which are UV attracted by the GFP but there are 
certain ones which are UV repulsed. In the latter case the perturbative flow equations give that the corresponding 
couplings go to infinity, which seems nonphysical. However the model can be safe from divergences if besides the GFP 
there exists a further nontrivial fixed point in the finite region of the phase space which UV attracts the trajectories, 
as the GFP does in the perturbatively renormalizable case. It can make the physical quantities finite, due to the 
finiteness of the dimensionless couplings. The nontrivial fixed point is usually referred to as a non Gaussian fixed 
point (NGFP). This is the main idea of asymptotic safety, Shortly, it means that there exists an UV attractive NGFP 
in the phase space. 

The paper is organized as follows. In Sect. II we give a short introduction on the functional RG method with some 
basic concepts and tools. Then the d-dimensional O(N) model is investigated as a classical example of the asymptotic 
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freedom. In Sect. Ill we treat the asymptotic safety, as the generalization of the idea of the asymptotic freedom. It is 
shown in the framework of the nonlinear a, Gross- Neveu, sine Gordon and quantum Einstein gravity models how the 
asymptotic safety takes place. In Sect. IV the conclusions are drawn up. 

II. RENORMALIZATION 

The functional renormalization group method is one of the most useful nonperturbative tool to investigate quantum 
field theoretic models. We define a model by its action at a high (UV) energy scale, because at UV the form of the 
action is generally quite simple due to the high symmetries appearing at large energies, which belongs to the short 
distance microscopic interactions. The UV action contains interaction terms multiplied by the UV couplings. The 
RG method systematically eliminates the degrees of freedom of the theory in order of decreasing energy, and provides 
us the value of the couplings at lower energy scales. As a result we obtain the value of IR couplings at practically 
zero energies. The IR the effective potential usually has an involved structure. There the long-range interactions may 
induce further couplings and can induce non-localities, or global condensates. 

The nonperturbative nature of the functional RG can be nicely enlightened, if one considers the tunneling, which is 
a highly nonperturbative phenomenon in Quantum Mechanics. In the path integral formulation of QFT the extremum 
of the UV action is the classical path. Thus if one treats QM then the tree level RG evolution corresponds to this 
classical path. It implies that the RG method should dress up the couplings in such a way which can drive the 
classical path to a path which corresponds to the tunneling process. Naturally it requires all the modes till the deep 
IR regime [8-12]. Moreover we need further terms in the gradient expansion. We note that the tunneling can be 
treated by a more exact manner in the framework of solving Schrodingcr equation with the corresponding potential 
numerically. However the Schrodinger equation can be considered as an effective theory, while the path integral 
formulation of QFT is more fundamental. The wavefunction which corresponds to the tunneling should be combined 
from practically plane waves since the field variables in QFT are expanded from them. If one includes higher order 
terms in the derivative expansion then more involved functions can appear in the field expansion. This fact requires 
to take into account all the modes and the wavefunction renormalization. Both can be systematically treated in the 
framework of the RG method. 

The quantized anharmonic oscillator coupled to a heat bath provides us a system, where the "holy grail" of the phase 
transitions, namely the quantum-classical transition can be investigated [13-18]. The heat bath can be considered 
as the environment of the original quantized oscillator, which can be identified by the system. The environment 
can be represented by harmonic oscillators which are coupled to the system. The RG procedure integrates out the 
environment and leads to an effective system. It appears as if the quantum effects are dissipated by the environment 
driving the system to a classical regime. 

The RG method can be used in a powerful manner in almost every area of modern physics. In condensed matter 
systems, e.g. one can investigate the Bose-Einstein condensate - BCS superconductor transition for ultracold fermionic 
atoms [19-21]. Furthermore the RG method can account for the essential scaling of the correlation length, which 
typically appears in low dimensional thin superfluid film structures. The appearing Kosterlitz-Thouless (KT) type or 
infinite order phase transitions [22, 23] can be described in the framework of the 2d 0(2) model [24, 25], in the 2d 
sine-Gordon (SG) model, too [26-29], or its generalization in fermionic models [30, 31]. 

New ideas or improvements of the RG method are usually investigated in the d-dimensional O(N) model. The 
calculation of the critical exponents in the 3d O(N) model provides us a widely acceptable testing ground of the 
functional renormalization. One can get the exponents by field expansion of the potential [32-34] , or by e-expansion 
in 4 — e dimensions [35-37] . The convergence of the exponents in the derivative expansion is also investigated [38-42] , 
however preciser results can be obtained without expanding the potential [43-45]. Furthermore, the supersymmetric 
version of the model has also been attracted a considerable attention [46, 47]. The exponents can be determined from 
IR scalings, too [48]. Likewise, less computational effort can give very good results for the exponents by using the 
BMW approximation, where the full momentum dependence of the correlation functions is considered [49-51]. 

The RG treatment of 2+1 fermionic systems has a revival since the new results on graphene. The Gross- Neveu 
model [52-54] give on the other hand a good testing ground to investigate a model with NGFP and asymptotic safety. 
The Thirring model [55, 56] serves as a good playing ground to investigate the chiral symmetry breaking, which 
appears in electroweak interactions. The nonperturbative treatment of the RG method is capable of describing quan- 
tum electrodynamics [57-59], non-abelian gauge theories [60-67], furthermore the confining mechanism in quantum 
chromodynamics even in finite temperature and chemical potential [68-71]. 

The nonlinear a model itself has a widebroad interest in many branches in quantum physics starting from phe- 
nomenological aspects of high energy physics, condensed matter systems, and strings. Furthermore it exhibits a 
symmetric and a broken symmetric phase, similarly to the O(N) model and the quantum Einstein gravity (QEG), 
and the model also has a NGFP in the UV [72-76] . The other common feature with QEG is that they have nonpoly- 
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nomial interactions, the couplings have the same dimension there, and both models needs background field technique 
to construct the RG equations [74] which is widely used in gauge theories [77] and in QEG [78]. Among the low 
dimensional scalar models we also mention the 2d SG model [79], which is non-trivial quantum field theory with 
compact variables - similarly to the non-Abelian gauge theories - which is supposed to be the key to the confinement 
mechanism. The functional RG treatment shows that the SG model has two phases [28, 29, 80, 81] and there is an 
IR fixed point in the broken phase [28, 29], furthermore the model has both an UV Gaussian and a NGFP, and what 
is more, the latter shows singularity. 

The IR limit of the RG flows is also a great challenge to reach. In scalar models the IR scaling of the symmetric 
phase can be easily obtained. It can be characterized by such an effective potential that has a well defined single 
minimum at the origin at the value of the field variable <j> — 0. In the phase with spontaneously broken symmetry 
the effective potential tends to be degenerate which means that a plateau starts to form in the potential around the 
origin [2, 6, 82]. It is due to the huge amount of soft (nearly zero energy) modes which can excite the ground state 
without energy practically [3, 83, 84]. In this region of the phase space the perturbative treatments naturally do not 
work. However the RG method uncovers us that there is an IR fixed point in the broken phase [28-30, 82, 85]. One 
can define the correlation length in the IR regime which can help us to determine the corresponding critical exponent 
v, moreover one can also determine the order of the phase transition of the model [48, 86, 87]. The IR fixed point 
explicitly shows the limitation of the theory where it can be treated by its original degrees of freedom. The bulk 
amount of soft modes show that new elementary excitations arises in the model at low energies and its treatment 
needs a new model, or at least some new interaction terms. 

The other limit towards the UV gives a further challenge of the RG method. Recently the UV limit of QEG is 
widely investigated [78, 88-93]. The model usually has a GFP. However around the GFP the Newton constant or 
coupling starts to blow up in the UV limit, which makes the physical quantities infinitely large there. It implies that 
that the model is perturbatively nonrenormalizable [94, 95]. As a possible solution for this problem it was conjectured 
[96-99] and later showed in low dimensions [100, 101], that there is a further UV NGFP fixed point in the phase space 
which makes the Newton constant finite and the model becomes safe from divergences, which is called asymptotic 
safety [91, 102-105]. Interestingly the model also has an IR fixed point in its broken phase [87, 106-110]. 



A. Blocking 

QFT can be formulated by using the path integral formalism [111] . The generating functional expresses the vacuum- 
vacuum transition amplitude and has the form 

Z = j V<pe~ Sk = J d<j) . . .d<j) k - A kd(f>kd(f)k+Ak ■ ■ ■d(f) 00 e~ Sk (1) 

where S k is the (blocked) action at the momentum scale k. The extremum of the action gives the classical path. The 
path integral is performed for all the possible paths between the given initial and final states. In order to handle them 
we order the paths according to e.g. by their momentum [112] or even by their amplitude in the internal space [113]. 
The procedure leads to the Wetterich equation [1] for the effective average action [114]. We can assign a certain k 
to each path, usually we enumerate them by their decreasing value. We note that the momentum k is basically a 
bookkeeping device for the modes, and it has no direct connection to any value of energy or momentum of the modes. 
The RG method provides us a systematic treatment to take into account the quantum fluctuation systematically. The 
modes are integrated out one by one in with decreasing value of k. After the elimination of some modes, the action 
changes to Sk-Ak, which is given at the lower scale k — Ak. If one eliminates the modes in the vicinity of the UV 
scale, then we recover the perturbative RG treatment. 

The highest scale k is the UV cutoff, which is denoted by k\. The blocked action is built up on e.g. symmetry 
considerations or analogies to other field theoretic models. Usually the action contains a kinetic and a potential 
terms, the latter contains interacting terms that are multiplied by couplings, and they carry the scale dependence in 
the action. The couplings usually have physical meaning, they can be related to e.g. particle masses or interaction 
couplings. The initial value of the couplings should be given at a high energy UV scale describing short-range 
interactions. If one integrates out the modes between the scales k and k — Ak then the value of the coupling changes. 
It is illustrated in Fig. 1, where the coupling is denoted by g k - 



The figure illustrates how the value of the coupling changes as the scale k decreases. The scale k starts at k\ and 
tends to zero. One has to integrate out infinitely many degrees of freedom, where the modes are enumerated by the 
continuous index k. The coupling becomes a scale dependent function, the RG evolution eventually gives how the 



(IR) <- k 



k — A.k k k + Ak 



k^oc (UV) 



FIG. 1: The change of the coupling g is presented as the scale changes from k + Ak — > k — > k — Afc. The limit k — > oc 
(k — ¥ 0) corresponds to the UV (IR) limits, respectively. 

value of the coupling changes till we reach the scale k = 0, and then we obtain the effective potential which contains 
the coupling go. From this point of view the functional RG method can be considered as a tool to perform the path 
integral in the generating functional. 

The RG method gives a partial differential equation for the action, however if one takes a functional ansatz for its 
functional form, e.g. 

Sk[4>,9i]=Y,9i(k)J r i(4>), (2) 

i 

then one can deduce a system of differential equations for the couplings. Here k denotes the scale, g%{= gi(k)) are 
the dimensionful couplings, and <ft is the field variable. The functionals T is typically Taylor expanded according to 
its field variable but there are situations where a Fourier expansion is required if one would like keep the periodic 
symmetry. 

B. The Wetterich equation 

The RG evolution equation for the blocked action is derived by Wilson [115] and Polchinski [116]. Wilson's idea 
considers the quantum fluctuation by decreasing order in the scale k, which can correspond to the summation of the 
loop expansion. There the short distance interactions are integrated out in order to get an effective long distance 
theory. Polchinski's RG method seems to sum up the perturbative expansion systematically [3]. 

The Wetterich equation is a functional integro-differential equation for the effective action [1, 2}. We start its 
derivation from the generating functional 

Z = e w ^ = J x»[0] e -( s *+ w *W- J -*), (3) 

in Euclidean spacetime, where J denotes the source of the field variable. We use the notation / • g = d d xf(x)g(x). 

The integral gives the average of fields over a volume k^ 1 . One integrates out the modes with the scale larger than 
k. There is an additional term besides the action, which is called the IR regulator, and it has the form 

n k [4>] = ^4>-n k -<(>, (4) 

which acts as an IR cutoff. Here the dots also denote momentum integrals. The IR regulator guarantees that the 
modes with its scale larger than k are taken into account in an unaltered way while the modes with scale smaller than 
k are suppressed. It should satisfy the conditions 

1. lim TZk > 0, i.e. it serves as an IR regulator, since it removes the IR divergences, 

p 2 /fc 2 ^0 

2. lim TZk — > 0, which expresses that if the regulator is removed, then we get back Z when k — > 0, so we obtain 

fe 2 /p 2 ^0 

back the original model in this limit, 

3. lim TZk ~ * 00 j which ensures that the microscopic action can be recovered in the limit S — lim^A IV It also 

k 2 — >oo 

serves as an UV regulator. 
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We have many possibilities to choose a functional from to the IR regulator [4, 34, 117-120]. Some typical regulators 
which are used frequently in the literature are 



n k 



-b{p 2 /k 2 



ap 



1 - e -6(p 2 A 2 ) c 



(exponential), 



TZk = p (power law), 

T^k = (k 2 — p 2 )0(k 2 — p 2 ) (optimized or Litim's). 



(5) 

(6) 
(7) 



The effective potential should not depend on the IR regulator but due to the truncation of the functional ansatz 
for the action induces regulator dependence, therefore one has to check how the obtained flows depend on them. 
Differentiating both sides of Eq. (3) by k one obtains 



d k W k [J] = e- w *^ J V[4>]d k K k e 



= -e- w ^d k TZ k 



_S_ 

SJ 



-(S k +K k [4>]-J-<t>) 



,W k [J] 



The effective action T k [<j>] is the Legendre transform of W k [J] , i.e. 

T k [</>] = -W k [J] + J- 

with the field variable 

5W[J] 



sj ■ 



The derivative of the effective action T k [(f>] w.r.t the scale k is 



dkT k [<f>] = -dkW k [J]- 



SW[J] 
SJ 



d k J + d k J<j>=-d k W k [J]. 



(8) 



(9) 



(10) 



(11) 



We redefine the effective action according to T k [(j>] + TZ k [4>] — > Tk, introduce the 'RG time' t = log ^ and then we 
obtain the Wetterich equation 



l Tr . K k 



2 K k + T'l 



(12) 



where the notations ' = d/d(p and ' = d/dt are used and the trace Tr stands for the integration over all momenta 
and the summation over the internal indices. The functional form of the effective action is assumed to be similar to 
the microscopic action 



= 5>(*ow). 



(13) 



Typically the effective action for scalar fields has the form 

1 



z k {<p x ){9M 2 + Wx) 



(14) 



where Z k = Z k (<f) x ) = Z k (<f>,p) is the wavefunction renormalization. It is the next to leading order contribution to 
the gradient expansion of the action in the field variable after the local potential approximation (LPA) . When we use 
the latter approximation the wavefunction renormalization does not evolve, implying that Z k = 1. Further terms in 
the derivative expansion can be 



d d x 



v k (<f> x ) + \z k (c» x )(d^ x ) 2 + ffiOfeXd^) 4 + H 2 (c» x )(n<p x ) 2 + ... 



(15) 
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If one inserts the form of the effective action in Eq. (14) into Eq. (12), then one obtains the following evolution equation 
for the potential 



V k = 



1 



2 J p z kP 2 + n k + v k " 



(16) 



with the d-dimensional momentum integral. Assuming that the wavefunction renormalization is momentum indepen- 
dent, i.e. Zk(4>,p) = Zk{4>) = Z k , then one obtains the evolution equation 



Zk 



^z' 2 p 2 + 4z' k (z' kP 2 + vn |pW + vn 2 (z k + d p ,n k ) 



-2 



[ P 2 z k + n k + v k "] 2 ( P 2 z k + n k + v k "f { P 2 z k + n k + v k "f 

{z' kV 2 + v>») 2 (z k + d p .n k + l P 2 dl 2 K k ) + \Z' k p 2 {z' kP 2 + v?) (z k + d p .n k ) 



(p 2 Z k +n k + v^y 



(17) 



for the wavefunction renormalization. 

If one Taylor expands the potential V k by its field variable, then one arrives at the potential for d-dimcnsional one 
component scalar (j> 4 model of the form 



N 



v k = j: 



92n ,2r. 



(18) 



After inserting it into Eq. (16) then one obtains a system of ordinary differential equations for the evolution of the 
couplings. By using the Litim's regulator in Eq. (7) the momentum integral in Eq. (16) can be analytically performed 
in any dimensions. It is 



V k = 2v d k' 



,2 k 



dk 2 + V" ' 



where 



Vd 



2 d + 1 n d / 2 T(d/2)' 

with T(d) the Gamma function. The form of the evolution equations for the dimensionful couplings is 

9i = Pi(9j,k), 

with the P functions. In case of the of the <j> 4 model their general forms are 

Ilk 



2 + n k + v" 



in LPA, which becomes 



Pi{ gj ,h)=dl [2v d k d - 



dk 2 + V 



(19) 



(20) 



(21) 



(22) 



(23) 



if one uses the Litim's regulator. 



C. Evolution equations 

One can reformulate the evolution equations in Eq. (21) and the (3 functions into dimensionless expressions according 

to 

9% = -di9i + ai{gj) = (24) 

where cti(gj) — fii{gjk~ d i ,1). The dimensionless couplings can be related to the dimensionful ones as <ji = k^ dj g^ 
with d the canonical (mass) dimension. 
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The evolution equations are usually highly nonlinear system of ordinary differential equations, which have no 
analytic solutions in general. If the phase space which is spanned by the dimcnsionless couplings is of high dimension, 
then it is extremely difficult map out the whole phase structure numerically. Therefore it can be useful to find such 
a tool which enables us not to avoid the most important parts of the phase space, i.e. where the evolution slows 
down or even stops. The points where the latter situation takes place are called fixed points. The fixed point of the 
evolution equations is defined as the singularity or stationary points of the evolution equations, i.e. which satisfies 

iji = 0. (25) 

The couplings which are the solutions of the system of algebraic equations in Eq. (25) are denoted by g* . Naturally 
the fixed points can be usually found only numerically. In the vicinity of the fixed points, due to the slowing down of 
the RG evolution, the flow equations can be linearized. We note that the linearized flow equations can differ around 
different fixed points. If we introduce yt = gt — g* we obtain the linearized evolution equations 

Vi = Mijyj, (26) 

with the matrix 

M l3 = ^, (27) 

that can be constructed by taking the derivative of the j3 functions w.r.t. the dimcnsionless couplings. The eigenvalues 
of the matrix M are denoted by s n . One can diagonalize M, with the help of a linear transformation represented 
by the matrix S which satisfies the relation S^ k M/~iSi n = 8 in s n . We introduce z% = S~^yk in order to decouple the 
linearized flow equations according to 

z t = s^. (28) 

Its solution reads as 

Zl = Zi(k A )e Btt = Zi(0) (^-) ' , (29) 



k A 

where k A is some reference scale, depending on the fixed point under investigation it can be e.g. the UV cutoff. The 
real part of the eigenvalue determines whether the trajectory is attracted or repelled by the corresponding fixed point. 



D. Classification of fixed points 



Let us assume that we have two dimensionless couplings, thus the phase space is 2-dimensional. It is quite straight- 
forward to generalize the classification for higher dimensional phase spaces, although the structure is a bit more 
complicated. It is clear that the signs of the eigenvalues determine whether we approach or go away from the fixed 
point where the linearization is performed. In this simple example case the eigenvalues are the solutions of a second 
order algebraic equation, therefore in general they are complex numbers. The imaginary part of the eigenvalues cannot 
alter the distance from the fixed point, altogether they make some oscillations around it. Thus the real parts of the 
eigenvalues determine the types of the scalings of the couplings. 

Let us denote the eigenvalues by s\ and s 2 , and consider the UV limits, i.e. k — > oo. According to the solution in 
Eq. (29) we have six possibilities. 

1. The eigenvalues are real, Si,s 2 € R and they are negative, Si,s 2 < 0. Then the trajectory approaches the 
fixed point, and is called an attractive fixed point. 

2. si, s 2 € M and s\, s 2 > 0. Then the trajectory goes away from the fixed point, and it is called a repulsive fixed 
point. 

3. si,S2 € M and with opposite signs. Then there is a direction which flows into the fixed point, and there is 
another one where the flow is repelled by the fixed point. The fixed point is called a hyperbolic point or a saddle 
point. 

4. The eigenvalues are complex Si,s 2 € C. They are necessarily constitute complex conjugate pairs. Here the 
signs of the real parts of the eigenvalues determine how the trajectories behave. If Sftsi, Ks 2 < 0, then the fixed 
point attracts the trajectories. The imaginary parts of the eigenvalues give some oscillation for the trajectory. 
The fixed point is called an attractive focal point. 
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FIG. 2: The irrelevant and relevant scalings of the couplings are shown. In the irrelevant case the the value of the coupling 
tends to oo with a power law behavior as k -¥ 00, and the scaling remains unchanged then, since the linearization gives better 
and better approximation. The value of the couplings decreases when the trajectory goes away from the UV regime. Its power 
law behavior is limited to the attractive region of the fixed point, far from its scaling regime it can deviate from the power 
law one. In this example it appears at k ~ 0.1. Then a new linear section can appear due to a possible other fixed point in 
the phase space. Naturally the new fixed point can turn the irrelevant scaling even into a relevant one. In the other figure the 
relevant scaling is shown. There the g — > as k — ► 00, while it grows up in the opposite direction of the scale. 



5. Si, S2 € C and z 5isi, 3?S2 > 0. Then the trajectory is repelled by the fixed point, and it is called a repulsive 
focal point. 

6. Si,S2 € C and the real part is zero. Then, due to the oscillation the trajectory circulates around the fixed 
point along a closed trajectory. It is called an elliptic point. 

Here we assumed that the eigenvalues are nonzero. If one of them is zero, then the other determines, whether the 
trajectory is attracted or repelled by the fixed point. If both are zero, then one should calculate those terms beyond 
the linearized approximation in the (3 functions in Eq. (24), which can give nontrivial contribution. We also note that 
he classification according to the IR limit k — > using Eq. (29) gives opposite trends for the trajectories. 



E. The Gaussian fixed point 

The GFP corresponds to the origin of the phase space, i.e. g* = 0. It describes a free theory for massless particles. 
Let us see the linearization of the flow equations in the vicinity of the GFP. If one Taylor expands the j3 functions 
around the origin of the phase space then one obtains that 

Pi = -difji + aigi + a lok g g k . . . (30) 

The matrix M defined in Eq. (27) is 

Mij = -dij + a,ij, (31) 

but one can prove that = 0, when i < j, therefore the eigenvalues are simply equal to the negative of the canonical 
dimensions, i.e. 

= -di. (32) 

It implies that all the eigenvalues are real in the GFP. The sign of the eigenvalue determines how the coupling behaves 
as we approach or go away from the GFP. If Sj > which means that the canonical dimension is negative di < 
then the transformed coupling Zi — > 00 when the scale k — > 00 or the RG time t — > 00, implying that gi — > 00 so the 
trajectory is repelled by the GFP. We illustrate it in Fig. 2. If the scale k is lowered then the coupling tends to zero, 
and it becomes less and less important as k — > 0, therefore we call the coupling gi as an irrelevant one. It implies 
that the irrelevancy is defined according to the IR scaling behavior, and it represents how the flow changes if the 
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trajectory goes away from the fixed point. Otherwise this definition is meaningful only for the GFP, since it assumes 
real eigenvalues. Its generalization is however quite straightforward by considering the real parts of the eigenvalues. 
One cannot conclude from an UV behavior to the IR one, since the trajectory can approach many other fixed points 
which can totally overwrite its scaling behavior there. 

The other possibility is when the eigenvalue of M in Eq. (27) is negative, i.e. Si < giving positive canonical 
dimension di > 0. Then Zi —> i.e. gi — > g* when k or t — >■ oo. In this situation the trajectory is attracted by the 
GFP, see Fig. 2. The coupling gi is said to be relevant, since it becomes more and more important as the trajectory 
goes away from the fixed point. The appearance of a new fixed point besides the GFP can overwrite the scaling 
properties of the couplings. Furthermore they can save the flows from the singularities [29, 48, 86, 87]. A coupling is 
said to be marginal if the corresponding eigenvalue is zero. The dimension of the critical surface equals the number 
of negative eigenvalues. 

We call a theory perturbatively renormalizable if V di > in the vicinity of the GFP. In this case we have only 
renormalizable couplings. A theory is asymptotically free if every value of the coupling tends to zero in the UV 
limit, i.e. lim gi — for V gi around the GFP. Examples for asymptotically free theories are e.g. QCD, and the 

t— >oo 

3-dimensional </> 4 model. 

The values of the irrelevant couplings blow up in the UV, and this makes measurable quantities nonphysical, since 
they will be infinitely large. However one can relate them to critical exponents. If s is the negative eigenvalue of the 
matrix M, then its negative reciprocal gives the mass critical exponent v, or the exponent of the correlation length 
£, i.e. 

v = -\js. (33) 
Let us recall that the scaling of the correlation length is 

S~(T-T C )- V , (34) 
for continuous or second order phase transitions, where T is the temperature, and T c is the critical temperature. 



F. d-dimensional O(N) model 

As an example of asymptotically free theory, we treat the d-dimensional O(N) model by functional RG method. We 
map out its phase structure, find the fixed points and the corresponding exponents. Although most scalar models are 
only toy models, but the O(N) model has experimental realization for certain values of N, i.e. they can characterize 
the following physical systems: 

N = polymers, 

N = 1 liquid-vapour transition, or uniaxial (Ising) ferromagnets, 

N = 2 He 2 superfliud phase transition, 

N = 3 Heisenberg ferromagnets, 

N = 4 chiral phase transition for two quark flavors. 

The 3d 0(1) or 3d </> 4 model is a widely investigated model, and it possesses a nontrivial fixed point. Using Eq. (19) 
for d = 3, one can get the following evolution equation for the dimensionful potential 

1 fc 5 

^^W' (35) 

if we use the Litim's regulator. If the functional form of the potential in Eq. (18) is inserted to the evolution equation, 
then one can obtain the /3 functions for the couplings. Their dimcnsionlcss forms are 

h = -2.92 - ~ 9A 



87r(l+.9 2 )V2' 

h = ~~ 9i+ 16n(l + g 2 )^ (36) 

where we neglected the evolution of further couplings which are generated by the RG procedure. In the symmetric 
phase their effects are really negligible, but they play crucial role in the evolution broken phase. Other regulators 
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would give qualitatively similar equations but e.g. with different multiplication factors, or different exponents in the 
denominator in the second term on the r.h.s.. Here we do not derive the perturbative RG equations, instead we 
expand the (3 functions in Eq. (36) by the coupling g 4 . The zeroth order approximation gives the flows driven by the 
canonical dimensions. They are 

g i =g i (k A )e- dlt , (37) 

with gi(k\) the initial value and 

di = d + i(l-^J , (38) 

the canonical dimension of the coupling. This approximation of the flow equations has a GFP. The matrix M in 
Eq. (27) is 

with the eigenvalues s\ = —2 and S2 = — 1. Thus the GFP is an UV attractive point, and both couplings are relevant. 
The eigenvalues and the type of the GFP do not change if one considers further terms in the approximation. Taking 
into account the next term in the expansion of the flow equations in Eq. (36) one gets 

.94 = -54 + H + 0{gf). (40) 

The structure of the phase space is plotted in Fig. 3. The model has two fixed points. They can be found by solving 
the system of algebraic equations fa = and fa = 0. The first is the GFP, g^ = g% = 0, with the eigenvalues s± = —2 
and S2 = — 1 found at the tree level previously. Furthermore, there is a nontrivial fixed point at g\ = —1/3 and 
g% = 167r/3. The matrix M now has the form 

M -("o 2 -r4) ss= -v3 W = ("o 2 ^) (4I) 

Its eigenvalues are si = —2 and S2 — 1, so the fixed point is a saddle point or a hyperbolic point. It is called 
Wilson-Fisher (WF) fixed point. It appears in 4> A models with dimension 2 < d < 4. When d — > 4 then the WF fixed 
point tends to the origin and in d = 4 it melts into the GFP. Since the critical exponent v of the correlation length £ 
is identified as the negative reciprocal of the single negative eigenvalue of the matrix M coming from the linearization 
of the evolution equations, then it gives v = —l/s\ = 1/2. The approximation of the model with two couplings makes 
the problem a mean field type one. 
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FIG. 4: The phase structure of the 3d cj> 4 model calculated from the flow equations in Eq. (36). The model has a GFP, a WF 
and an IR fixed point. The latter one corresponds to a crossover fixed point between the UV and the IR ones. 



The 3d (f> 4 model has two phases. The trajectories tending right from the WF fixed point in Fig. 3 correspond to 
the symmetric phase. The trajectories tending left belong to the broken phase. There the Z 2 symmetry (<j> <H- —(f)) is 
broken. 

Further terms in the expansion in g^ do not give further fixed points, but only changes qualitatively the position of 
the WF fixed point. If one considers the exact RG flow equations in the two couplings in Eq. (36) we have the GFP 
at the origin and the WF fixed point at g\ = — 1/4 and g% = 2v / 3~7r. The matrix M is 



Its eigenvalues are s\ = — 2 and s 2 = 4/3. The WF fixed point remains a saddle point. The higher order terms in 
approximations (e.g. consideration of further couplings, inclusion wavefunction renormalization or even higher order 
terms in the derivative expansion) cannot change the type of the fixed point only the position is shifted further. We 
illustrate the phase structure in Fig. 4. Interestingly the trajectories tend to a single point in the broken phase. In 
might assume that there can be a further fixed point in the IR limit. That seems strange according to the j3 functions 
in Eq. (36), since we can find all the fixed point solutions analytically. The only possibility is that the ratios in Eq. (36) 
can tend to zero in the way that both the numerator and the denominator vanishes. If they are of the form of 0/0, 
then it implies that the values of the couplings are g\ = — 1 and g\ =0. It corresponds to an effective potential of 
the form Vq = —(f> 2 /2. Clearly this fixed point cannot be reached by any expansions around the GFP since this point 
is a singular point of the evolution equations. 

The effective potential of the symmetric phase can be handled quite easily numerically in contrast to the one in 
the broken phase, where one can face a very difficult numerical and theoretic problem [2]. The reason is that there 
is a condensate in the broken phase [81, 84, 121-124]. It can be considered as a macroscopic object constituting 
of a huge amount of soft modes. Its measure gives a dynamically appearing momentum scale in the model, which 
implies that there cannot be such modes which momentum is that is larger than the one characterizing the condensate 
[84, 114, 125]. The evolution stops when one reaches the scale of the condensate, which manifests in the form that the 
flows arrive at the singularity. Is there really a singularity there, or it is only a numerical artifact? It is known that 
the Wetterich equation in Eq. (12) should not drive the flows into singularity. However the approximation due to the 
strong truncation of the functional ansatz for the effective action, or in the derivative expansion may induce a singular 
behavior suggesting that the singularity is a numerical artifact. Therefore it is argued that the singularity can be 
avoided by a proper choice of the IR regulator, and the flow can reach arbitrarily small scales. However it was shown 
for a huge set of scalar models [29, 48, 86, 87] that such a singularity possesses a specific scaling behavior induced 
by the IR fixed point, therefore it may have a significant physical importance. The singularity in the RG evolution 
is always reached, but the RG equations do not loose their validity, they simply stop at a finite scale [44, 126]. This 
finite scale also appears when it seems that the evolution avoids the singularity with special regulators. In those 
cases, after an abrupt change in the scaling of the couplings, a marginal scaling appears towards the deep IR regime 
without any universality. Therefore the dynamically appearing finite momentum scale seems to be the only universal 
behavior in the IR. 




(42) 
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The dynamical scale is induced by an IR fixed point in the broken phase. One can determine numerically the 
critical exponents in its vicinity. The IR fixed point can be uncovered analytically from the /3 function in Eq. (36) by 
reparametrization of the couplings according to co = 1 + g 2} x = .94 /"^ an d d T = uid t - We obtain that 

d r Lo — 2uj(1 — u) — 

drX = ~X+£- (43) 

The reparametrized flow equations give the Gaussian (wj = 1, \*g = 0); an d the WF (wJ^ F = 3/4, XiVf = ^ 7r ) nxe d 
points, however another one appears at to* IR — and X*ir = 4tt. The latter can be identified with the IR fixed point 
where the trajectories of the broken phase meet. The corresponding eigenvalues are s\ — 1 and s 2 = 3/2 expressing 
the UV repulsive (IR attractive) nature of the IR fixed point. 

The determination of the critical exponents in the vicinity of the IR fixed point is demonstrated in the framework 
of the d-dimcnsional O(N) model. The potential has the form 

Nx A 



(p-K) n , (44) 

n=2 



n! 



with N\ the degree of the Taylor expansion and the dimensionless couplings k and A n for n > 2. For shorthand we 
use A 2 = A. The further dimensionless quantities are denoted by ~. The introduction of the coupling k serves a 
better convergence in the broken symmetric phase. The evolution equation for the potential reads as 



; / 1'-,/ (\ rj \ ( 1 , N ■ 



kd k V = -dV + (d-2 + r\)pV' + — r 1 - -rh; = ~~ + =~ , (45) 

1 d \ d + 2j \l + V' + 2~pV" 1 + V'J K ' 

by using Litim's regulator, with the notation ' = S/Sp. In Eq. (45) we introduced the anomalous dimension rj which 
is defined via the wavefunction renormalization according to r) = —d\ogZ/d\ogk and can be calculated by means of 
the couplings as 

" = — T+^a- (46) 

The inclusion of 77 in the RG equation mimics the evolution of the wavefunction renormalization. We note that a more 
precise treatment can be obtained if one Taylor expand the evolution equation for Z in Eq. (17). Again, we start with 
the 3d O(l) model. It is instructive to retreat the phase structure, since the IR fixed point moves to infinity. Now its 
flow equations are 



1 



K = —K 



27t 2 (1 + 2kA) 2 ' 



3A 2 

A = - A + ^TTW' (47) 

for the first two couplings if we set 77 = and A„ = for n > 2. The O(N) model in d = 3 has two phases. The 
typical phase structure is depicted in Fig. 5 for the couplings in Eq. (47). From the flow equations in Eq. (47) one can 
find two fixed points in the model. The Gaussian fixed point can be found at k* g = 1/2-k 2 and \* G = 0, although it is 
not situated in the origin due to the redefinition of the field variable. The linearization of the flow equations around 
the GFP gives a matrix with negative eigenvalues (sgi = — 1 and sq2 = — 1), i.e. the fixed point is repulsive or UV 
attractive. The WF fixed point can be found at the values of the couplings n* WF = 2/9tt 2 and X^p = 9ir 2 /8, with 
eigenvalues swfi = 1/3 and swf2 = —2, so it is a hyperbolic point. 

We usually identify the critical exponent v of the correlation length £ in the vicinity of the WF fixed point by taking 
the negative reciprocal of the single negative eigenvalue, which gives tjwf = 1/2 in this case. Let us notice that the 
flows tend to a single curve beyond the WF fixed point in the broken phase into the IR fixed point at k* ir -4- —00 
and A}^ = 0. Again, with the help of a rescaling of the couplings the attractive IR fixed point can be uncovered and 
one finds the following pair of evolution equations 

8 t lu = 2w(l-w)-^(3-4u;) ) 
U 2 

d T i = *(5w-6) + -5-(l-w), (48) 
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FIG. 5: The phase space of the 3d O(N) model is shown. The flows belonging to the symmetric (broken symmetric) phase 
tend to right (left), respectively. Again, the Wilson-Fisher fixed point plays the role of the crossover fixed point between the 
UV (which is now the Gaussian) and the IR ones. 




FIG. 6: The evolution of the anomalous dimension r\ is presented. In the left figure the flow of r\ shows a strong singularity 
as the function of the scale k, while this singularity becomes a power law like convergence as the function of the shifted scale 
k — k c in the right one. 



where u> = 1 + 2k\, £ = \/uj 3 and d T — kd^juj. The static equations now have the Gaussian (t* G — 0, luq — 1), the 
WF (£wf = 7I " 2 /3, ^wf = 3/2) and the IR (£* IR = 2ir 2 /3, u* IR = 0) fixed point solutions. Naturally the Gaussian 
and the WF ones has the same behavior as was obtained from direct calculations. However the new IR fixed point 
indeed corresponds to the values n* IR — > — oo and X* IR — 0, and the linearization in its vicinity gives the eigenvalues 
s 'im — ^ an d s' IR2 = 0, a positive and a zero one, showing that the fixed point is IR attractive, in accordance with 
the flows in Fig. 5. 

If we let many couplings and r] evolve, then the phase space does not change significantly, but the IR fixed point 
can be observed more easily. In order to demonstrate it we plotted the flow of the anomalous dimension in Fig. 6. We 
also plotted r\ as the function of the shifted scale k — k c in the right figure of Fig. 6. In the UV region the Gaussian 
fixed point drives the evolution of the anomalous dimension. In this regime it grows according to the power law 
scaling rj UV ~ (k — k c )~ 2 . There is a crossover (CO) scaling region between 10~ 8 < k — k c < 10~ 4 where a plateau 
appears giving a constant value for r\co ~ 0.043 due to the WF fixed point. Going further in the evolution towards 
the smaller values of k below k — k c ~ 10~ 8 one can find a third scaling regime. It appears a simple singularity in 
the left figure, but the shifted scale k — k c clearly uncovers the power law scaling of the anomalous dimension there 
according to r/m ~ {k — kc) 1 . This scaling region is induced by the IR fixed point. 

The evolution of the other couplings also shows such type of scaling regimes with similar singularity structure in 
the IR limit. There the power law scaling behaviors also take place as the function of k — k c with the corresponding 
exponents. This suggests that the appearing singularities are not artifacts and the RG flows can be traced till the 
value of k c . 

We note that one can find such a value of N\ where the evolution does not stop as in the upper figure in Fig. 6, but 
after a sharp fall during the flow of r/ it continues its RG evolution marginally giving a tiny value there. However, the 
singular-like fall possesses the same power-law like behavior as the function of k — k c , for any value of N\. It suggests 
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FIG. 7: The scaling of the correlation length £ as the function of the reduced temperature t, for various values of N. 
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TABLE I: The critical exponent v in the O(N) model for various values of N. 

that the value of rj rapidly falls to zero at k c and it is due to the numerical inaccuracy, whether the RG evolution 
survives the falling and can be traced to any value of A:, or the flows stop due to the appearing singularity. It strongly 
suggests that the singular behavior with its uncovered IR scaling for the shifted scale k — k c is not an artifact but is 
of physical importance, since the value of the scale k ci which appears there can characterize the condensate as the 
smallest available scale. If one tunes the initial values of the couplings Kk A , \k A etc., to the separatrix then the value 
of k c decreases. It enables us to define the correlation length £ in the IR scaling regime as the reciprocal of the scale 
where the evolution stops, i.e. £ = l/k c . As the initial couplings approach the separatrix the stopping scale k c — > 
therefore the correlation length diverges. Naturally it is infinite in the symmetric phase. We can identify the reduced 
temperature t in the O(N) model as the deviation of the UV coupling Hk A to its UV critical value, i.e. t <~ k£ a — Hk A - 
By starting evolutions for different values of the UV coupling Kk A then one can get the corresponding values £. One 
should fine tune the initial coupling to get larger and larger values of £ to reach the IR scaling regime. The critical 
initial value of k£ a can be determined by the trick, where one should fine tune its value in the log-log plot of the 
t, £ plane till one obtains a straight line there. This holds for continuous phase transition. The absolute value of the 
negative slope of the line provides us the exponent v corresponding to the correlation length. 

In the 3d O(N) model a second order or continuous phase transition appears according to the power law scaling 

£~t- v . (49) 

The results are plotted in Fig. 7. For a given value of N we obtain power- law like behavior for the scaling of £, and the 
slope of the line provides the exponent v in the log-log scale. The obtained results are listed in Table I. We denoted 
the WF (IR) values of v as (^ipj, respectively. The results show high coincidence. The values are taken 

from results obtained from derivative expansion up to the second order, since this approximation is the closest to our 
treatment. One can conclude that the exponent v can be also determined from the scaling around the IR fixed point, 
and has the same value as was obtained around the WF fixed point. 

In arbitrary dimensions one can consider the evolution of the couplings k and A around the GFP. It is situated at 
k* g = 2 1 - i TT- d / 2 {N + 2)/d(d - 2)T(d/2) and = 0. The corresponding matrix M is 

{ n j 3(W+2)4 1 - d 7r-" \ 

M = d-2r(i+d/2)2 ( 50 ) 

V d-4 J 

therefore the eigenvalues are sqi = 2 — d and sq2 = d — 4. When 2 < d < 4, then the GFP is UV attractive and 
the model is asymptotically free. However one of the eigenvalues become zero at d = 4, and the coupling A does not 



15 



necessarily converge to zero, and the model is not asymptotically free anymore. The appearing divergence can be 
related to the Landau pole, where a singularity of the couplings appears at a finite scale k, as e.g. in QED [58, 127]. 



III. ASYMPTOTIC SAFETY 



The concept of perturbative renormalizability and the asymptotic freedom is restricted to the GFP. It guarantees 
that the physical quantities which are calculated from the model do not suffer from divergences. However if there 
arc irrelevant couplings in the GFP which are crucial in describing the model that is investigated, (e.g. the Newton 
constant in QEG, which makes the GFP a hyperbolic one) then the systematic removal of the divergences may 
even induce infinitely many new important couplings which implies that the model disables us to give any physical 
predictions. 

However there can be further nontrivial (non-Gaussian) fixed points (NGFP), where the physically important 
couplings are relevant, which means that the theory can give finite physical quantities. This is the basic idea of the 
asymptotic safety. In general sense the asymptotic safety means that the theory is free from divergences if the cutoff is 
removed to infinity assuming that the corresponding fixed point possesses a finite number of UV attractive directions. 
In the fixed point those couplings should scale in relevant manner which are crucial to obtain finite physical quantities. 
Naturally there can be certain couplings which are generated by the RG procedure, but has no physical meaning. It 
is necessary to have finite number of relevant couplings, otherwise every trajectory would tend to the fixed point and 
the theory would not be predictive. The definition of the asymptotic safety is not restricted to the GFP. It implies 
that there is a UV NGFP in the phase space. The asymptotic safety requires that the eigenvalues corresponding to 
the linearized RG flows of the physically important couplings around the UV NGFP should have negative real parts. 

Examples for asymptotically safe theories are e.g. quantum gravity, the Gross-Neveu and the nonlinear a models 
with dimension 2 < d < 4, or the 2d sine-Gordon model. 



A. The nonlinear a model 



The nonlinear a model (NLSM) in general describes the dynamics of a map ip from a d-dimensional manifold M 
to a N-dimensional manifold Af. The model is renormalizable in d = 2 and is asymptotically free [128, 129]. However 
if one goes beyond d = 2 then it becomes nonrenormalizable since the existing UV attractive Gaussian fixed point 
becomes a hyperbolic one. However a nontrivial UV fixed point arises [128], which saves the UV limit, and the model 
becomes asymptotically safe. The action of the model contains only derivative interactions, e.g. 

S=\cJ d d xd^ a d»^h aP (v) (51) 

where h a p is the dimcnsionless metric, £ = 1/g^ and its dimension is [g ] = k^ d ^ 2 . The RG equations can be 
derived by background field method. The perturbative RG flow gives [72] 



d- 2_ R 

P§o = 

with 



Pgo = ^-^90 ~ c d ^9o, (52) 



Cd (4^)<v 2 r(d/2 + i)- (53) 

There are two fixed point of the flow equation in Eq. (52). At g^ G — we have a GFP, with the eigenvalue sg = 
(d — 2)/2. It is clearly positive for d > 2 giving irrelevant scalings and divergences in the UV limit. The other fixed 
point is situated at gQjy = (d — 2)N/ {2aN) giving the eigenvalue suv — 2 — d, which induces relevant scaling, 
therefore the NLSM exhibits asymptotic safety, yet at perturbative level. The critical exponent v of the correlation 
length £ equals the negative reciprocal of suv giving v = l/(d — 2). 
The functional RG flow equation results in the equation [72] 

P- = ^—^.9o - (54) 

2 1 - 2cd N(d+2)9o 

Again, we have two fixed points. The origin is a GFP with the same eigenvalue sg — (d — 2)/2, and we have a UV 
NGFP at g$j V = N(d 2 - 4)/(4c d dR). It is an UV attractive point with eigenvalue s uv = -2d(d- 2)/{d + 2). We 
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FIG. 8: The phase space of the nonlinear a model is presented. The trajectories tending to the left towards the IR fixed point 
are in the broken phase, the others belong to the symmetric phase. The fixed points are denoted by gray points, they are: UV 
attractive NGFP — > UV, hyperbolic NGFP — s- NG, and IR fixed point — > IR. The separatrix is situated between the UV and 
the NG fixed points, and is denoted by the thick line. The latter is a saddle point and plays the role of the crossover fixed 
point in this model. The dashed line shows the singularity limit, when go = —gi- 



note that this value is smaller then the one that was obtained perturbatively, therefore the exponent v is greater, 
e.g. it is v = 5/6, when d = 3. In the NLSM the perturbative and the exact RG flow equations give qualitatively 
similar results, and the functional RG method changes only the position of the fixed points, and the concrete value 
of exponents. The asymptotic safety appears in both the perturbative and the exact RG treatments. 

Taking into account more interacting terms in the action in Eq. (51) we should introduce further couplings in the 
model [130]. The /3 functions for two couplings reads as 

/%, = 9o+9o(N -2)Q d/2 . 2 + dg 1 (N -2)Q d/2+h2 , 
hi = -9i +9i(N-2)Q d/2t2 , 

where 

6 = 1 / (2n + 2 + d t )ff 20 + 4 + d t )g! \ 

( 47r )<V2r(n) U(« + !)(5o + ~9i) L n(n + 2)(g + ) 

The phase structure is given in Fig. 8. The model has two phases henceforward, but new fixed points appear. At 
9ong — 2/57T 2 and giNG = a hyperbolic NGFP appears. The trajectories are UV attracted from the direction of go 
and UV repelled from g\ . It implies that including a new coupling into the model the asymptotic safety obtained for 
a single coupling disappears. The eigenvalues are sjvgo = —6/5 and s^d — 2. The former gives v — — I/s^go = 5/6 
for the critical exponent, which coincides with the one obtained for g alone. However a further NGFP appears 
at g uvo = 16/357T 2 and gu V1 = — 12/357T 2 . The new fixed point is UV attractive as the flows show in Fig. 8, the 
eigenvalues are suvo = —0.457 and sjjvi — —13.11. Let us mention, that a further hyperbolic fixed point appears in 
the phase space far from the region shown in Fig. 8. Noticeably the trajectories meet at the GFP in the broken phase, 
if one follows the evolutions into the deep IR region, see Fig. 8. 

According to the (3 functions in Eq. (55) the quantity go + gi =0 makes the flows singular. We plotted it for 
different UV initial values of the couplings in Fig. 9. The trajectories have three different scaling regimes. In the UV 
there is a short, relevant scaling region induced by the UV attractive NGFP. Then the scalings become marginal in 
the vicinity of the hyperbolic NGFP. In the IR region the trajectories belonging to the symmetric phase diverge as k 1 . 
In the broken phase the flows tends to zero sharply at a certain scale k c , which can be identified with the correlation 
length £ again. Therefore one can determine the exponent v of £ in the IR. It gives v — 0.835, which is very close to 
the value v = 5/6, that was got analytically in the vicinity of the hyperbolic NGFP. 

B. The Gross-Neveu model 

The Gross-Neveu (GN) model is a 2-dimensional quantum field theoretic model of Nf flavors of massless relativistic 
fermions which interact via a four-fermionic term [131]. The model is asymptotically free. The GN model was also 
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FIG. 9: The flow of the quantity go + gi is shown. The full (dashed) line corresponds to evolutions in the broken (symmetric) 
phase, respectively. 



investigated at finite chemical potential in d = 2, where an analytic solution of the crystal ground state is found [132]. 
In d = 3 the model is not asymptotically free [85]. It can also be studied at finite temperature [52] and chemical 
potential [133]. The Euclidean effective action of the GN model has the form [85] 



S$,il>] 



(57) 



with the wavefunction renormalization which is set to Z$ = 1 in LPA. The dimensionless four-fermion coupling g 
can be obtained from the dimensionful g as 



g = Z-V-ig. 

The functional RG treatment of the GN model leads to the RG flow equation 

P g = (d - 2 + 2^)g - M 7 v d l[ (0)g 2 



(58) 



(59) 



in the Nf — > oo limit [85]. The constant d 1 denotes the dimension of the Dirac-algebra, e.g. it is d 7 = 4 when d = 3. 
The other constant Zf (0) is an IR regulator dependent quantity, for example Zf(0) = 2/cZ when one uses the Litim's 
regulator. The model has two fixed points. We have a GFP at g* = with the eigenvalue sq = d — 2. It shows that 
the four-fermion coupling multiplies an irrelevant operator, thus the model is perturbatively nonrenormalizable. The 
other NGFP is situated at 



.9 = 



4d. ( v d l[(0) ' 



(60) 



which, by using Litim's regulator in LPA and considering d — 3 becomes g* = 37r 2 /4. Its corresponding eigenvalue is 
suv = 2 — d, which gives a relevant coupling when d > 2, so this fixed point is a UV attractive NGFP. In the case of 
d = 2 the model is asymptotically free, perturbatively renormalizable with its UV attractive GFP. 

In the partially bosonized version of the GN model [52, 85, 134] a hyperbolic fixed point, and a UV attractive 
GFP appears, so the model becomes asymptotically free. From [85] the corresponding evolution equations for the 
d-dimcnsional bosonized GN model is 



-du +(d-2 + 7] a )u'p - 2d 1 v d l ( f )d {2h 2 p- + -^-2v d l*{u' + 2pu"- Vtr ), 



where we impose a polynomial ansatz for the potential 



oo » 



(61) 



(62) 



n=0 



18 



with the bosonic couplings A 2 „. One can relate the coupling g in Eq. (59) as g = h 2 /X 2 . The flow equation for the 
Yukawa coupling is 

h 2 = (d - 4 + 277V, + Va)h 2 + -^-Bvdh'li^iO, X 2 ;r hl> , r? CT ). (63) 



The anomalous dimensions are 



V* = S^-h z m\ (0;jfy), 

»7V = TT 8 ^ h2m h2 B)d (°> A 2;%,?7 CT ). (64) 
A// d 



We introduced the threshold functions 



d V d + 2/ 1 + w 

AFB)d( 2 1 // ^ \ 1 / f] a \ 1 

™4 (w;»^) - 7T— u + 7= ; ; ,,3 - o77~T + T ' 



(1+w) 4 d-2 (1 + u;) 3 V 2d - 4 4/ (1 + ^) 2 ' 

mi™^,^;^.^) = (^m) (1 + 0,0(1 + ^ - (65) 

Keeping only the couplings ft- 2 and A 2 we get the evolution equations for d = 3 



The flow equations have a UV attractive GFP at the origin at Kq — and X 2 q = and a non-Gaussian saddle point 
at hj^ G = 5.764 and X 2NG = 0.758 for iVy = 12. The phase structure can be seen in Fig. 10. The trajectories tending 
to left and right correspond to regions of different phases in the model. After the reparametrization of the couplings 




FIG. 10: The phase space of the 3d Gross-Neveu model is shown. The trajectories tending to the IR fixed point are in the 
broken phase, the others belong to the symmetric phase. The fixed points are denoted by gray points. The separatrix between 
the phases is represented by the thick line and connects the Gaussian (G) and the non-Gaussian (NG) fixed points. The GFP 
is a UV attractive, the hyperbolic NGFP is a crossover fixed point. 
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N 


1 


2 


6 


12 


50 


"IR 


0.922 


0.976 


0.990 


0.996 


1.00 



TABLE II: The critical exponent v in the 3d Gross-Neveu model for various values of N. 



according to uj = 1 + A 2 , x = h 2 /uj and d T = ojd t we get 

2 

d T uj = 2w(l-w) + ^-(5w-l), 

«* = 2x(*-2) + i ^(7 W + l-^). (67) 

The new flow equations have three physical fixed points. The UV attractive GFP can be found at Wq = 1 and = 0. 
The NGFP is situated at lo* ng — 1.758 and X*ng = 3-278 with a positive and a negative eigenvalues. There is a third 
fixed point at co} R = and \*ir — 355.206 with two positive eigenvalues implying that that this fixed point is a UV 
repulsive or IR attractive, therefore we can identify it with the IR fixed point at h? IR — and X^ir = —1 if we express 
them in terms of the original couplings. One can investigate the scaling of the couplings in the vicinity of the IR fixed 
point, and it is shown in Fig. 11. 
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FIG. 11: The scaling of the couplings in the 3d Gross-Neveu model is shown, x denotes hi and A2, the fixed points values are 
h*2 = and A2 = —1. The couplings have a long marginal scaling region induced by the hyperbolic NGFP. In the IR limit there 
is another scaling regime, which can be uncovered by using the shifted scale k — k c . 



The couplings are constants in the beginning of the flows, giving marginal scaling due to the hyperbolic NGFP in 
Fig. 10, that is x — x* <~ (k — k c )°. In the deep IR region the couplings have singular behaviors, i.e. they abruptly 
become zero or infinite at a certain scale k c . If one plots the couplings as the function of the shifted scale k — k c , 
then one can obtain a power-law like behavior according to x — x* <~ (k — kc) 1 , where x* is the fixed point value of 
the couplings, for example h\ = and A 2 = —1, see Fig. 11. It implies that the original four-fermion coupling g tends 
to zero in the IR limit. We also calculated the critical exponent v of the correlation length in the IR regime by the 
previously shown method, which is based on identifying the correlation length by the reciprocal of the stopping scale 
k c . The results are summarized in Table II. They show good coincidence with the ones obtained in [53]. 



C. The 2d sine-Gordon model 



The 2d sine-Gordon model can be an example for asymptotic freedom and asymptotic safety, simultaneously. Its 
effective action contains a sinusoidal potential of the form 



-(<9 M </>) 2 +u cos< 



(68) 



20 



where z is the field independent wavefunction renormalization and u is the coupling. The perturbative RG results 
beyond LPA [26] can account for the KT phase transition and give 1/z — >■ for 1/z < 8ir in the IR limit. By using 
the flow equation approach [135, 136] a different IR limit is obtained for the frequency, i.e. 1/z — > 4ir. However we 
note, that this method cannot recover the leading order perturbative UV results for 1/z < Air, due to the opposite 
sign obtained for the evolution equation for the frequency. Functional RG approaches can map the phase structure 
of the SG model in LPA [80, 81, 126, 137-139] and can also take into account the evolution of the wavefunction 
renormalization [28, 29, 140-142]. 

Besides the Z 2 symmetry the model has a periodic symmetry. The model has two phases. The effective potential 
should be convex (nonconcave), furthermore the RG equations should not break the periodic symmetry, too. This 
two requirements can be satisfied if the effective potential is flat, or zero [80, 143]. Then how one can distinguish the 
phases in the model? The simple answer is that one should consider the dimensionless effective potential, since only 
the dimensionful one should be convex. In the symmetric phase the dimensionless effective potential is flat, while it 
is a concave function in the broken phase. In LPA we got that it is Vo — —<fi 2 /2, which is repeated periodically in 
the field variable. This result can be obtained if one considers the evolution of higher harmonics in the potential, or 
if one follows the flow of the potential without its Fourier expansion. 

The phases cannot be always distinguished by the form of the effective potential. The real difference comes from 
the property how the effective potential depends on the UV initial value of the couplings. It seems to contradict to 
the fact that the effective potential is a simple function at k = 0. However this scale cannot be reached, since during 
the evolution we use that the ratio Ak/k is small disabling us to reach the exact value of k = 0. If we consider the RG 
time t then it is more apparent that we have only tendencies to infinity. Therefore it is meaningful to investigate how 
the asymptotic behavior of the effective action becomes the effective potential when the scale k is lowered. One can 
observe that a power law scaling appears in the deep region IR which can imply that the appearing scaling behavior 
can be extrapolated to k = 0. In the case of the 2d SG model the effective potential in the symmetric phase depends 
on the initial value of u(k\), while it is universal (independent on the initial coupling) in the broken phase. This type 
of separation of phases always works in any models. 



1. Local potential approximation 

By using the power law IR regulator in Eq. (6) with the choice 6 = 1 the leading order approximation of the RG 
evolution equations is 

u = —2u, 

i = 0. (69) 

These equation are usually referred to as UV scalings. The fixed points are u* — with arbitrary z*. These points 
constitute a line of fixed points. Since z does not evolve, then we are in LPA. The matrix M corresponding to the 
linearized RG equations in the vicinity of the line of fixed points is 

|~ 2 + ib ~4^\ = | 2 +i7^ °| (70) 

V >w V V 

The eigenvalue to the evolution of u is s = — 2 + j^r- Its sign depends on the value of z* . If z* < l/8w then s > 
and the coupling is irrelevant. When z* > l/8n then s < and the coupling scales in relevant manner. These two 
types of scalings gives the two phases of the model separated by the critical value z* = 1/8-k. The point u* = and 
z* = 1/8-k is called the Coleman point. The separatrix is a vertical line and goes through the Coleman point in the 
phase space. 

The symmetric phase, when z* < l/8n the phase is usually called the nonrenormalizable (or massless) phase due 
to the irrelevant scaling of the coupling u. The other phase is the broken one. There the Zi symmetry breaks down, 
it is called the renormalizable (or massive) phase (it is relevant there). 
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FIG. 12: The phase structure of the 2d sine-Gordon model treated in linear approximation. The separatrix is denoted by a 
thick line. 



2. Linearized RG equations 



By using the power law IR regulator in Eq. (6) with the choice 6=1 the linearized RG evolution equations are 

1 



The fixed point solution is u* — and z* arbitrary. The matrix M is 



-2 + — 



12tt 








(71) 



(72) 



In case of the linear approximation the wavefunction renormalization z does not evolve, similarly to the LPA case, 
and we also have the same eigenvalues. It also implies that in the vicinity of the line of fixed points z does not evolve, 
which is reflected by perpendicular trajectories to the horizontal axis. The eigenvalues are the same as was got in 
LPA. The phase space is presented in Fig. 12. The point z* = l/8n and u* c = is now called the Kosterlitz-Thouless 
(KT) fixed point. There is a separatrix which divides the phase space into two parts. The trajectories tending to the 
line of fixed points belong to the symmetric phase, all other trajectories constitute the broken phase. Around the KT 
point an essential scaling appears due to the infinite order phase transition in the model. 



3. Exact RG equations 



The RG equations can be obtained exactly when b = 1 in this approximation (where the potential contains only 
the fundamental mode, and the wavefunction renormalization is field independent). The flow equations are 



-2u + 



1 



2ttuz 
u 2 



24.7T (1 



£2)3/2' 



(73) 



The phase structure is plotted in Fig. 13. One can easily see that the phase structure of the linearized treatment 
in Fig. 12 is a part of this phase space around the KT point. The separatrix was a straight line there which is now 
curved. If one tries to find the fixed points of the model then immediately realizes that it has not got any. However 
the phase structure suggests, that there is the line of fixed points at u* = 0, and it seems that we have an IR fixed 
point in the broken phase at u* = 1 and 1/z* = 0. Furthermore a UV attractive NGFP appears at u* = 1 and z* = 0. 
Let us see them one by one. 



Line of fixed points 
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The line of fixed points can be found at u* = with arbitrary z*. Its scaling behavior can be obtained by expanding 
the RG flow equations in Eq. (73) in u which gives back the approximate equations in Eq. (71). Thus we have the 
same scaling behavior, i.e. a relevant scaling for z* > l/8n and an irrelevant one for z* < l/8n and the different 
scaling regimes are separated by the KT point at u* = and z* = 1/8tt. 

IR fixed point 

The IR fixed point is situated at u* = 1 and 1/z* — in the broken phase. It cannot be found directly from equations 
in Eq. (73), too, therefore we should reparameterize the couplings. After introducing lj = \/l — u 2 , \ = 1/zoj and 
d T = oj 2 kdk we arrive at the evolution equations 

2 

d T cu = 2lo(1-lu 2 ) - 

9tX = X2i 24^ " 2X(1 ~ " 2) + ~ W) ' (?4) 

possessing the fixed point w* = and arbitrary x* which corresponds to u* = and arbitrary z* , thus we get back 
the lines of fixed points. However there is another fixed point at \* — an d u * = Oj which can be identified by the 
IR fixed point at 1/z* — Q and u* = 1, which is IR attractive. 

If we introduce fc = min(zp 2 + R), the RG evolution becomes singular at k — k c when 

fc 2 -V fe "(</» = 0)| fe=fec =0, (75) 

where P = bk 2 [z/(b — l)] 1-1 / 6 , when 6=1, then k = k. The solution of this equation defines the scale where the 
potential becomes degenerate. 

The normalized coupling u tends to 1 for every value of b. It shows that the degenerate potential (which satisfies 
Eq. (75)) occurs in the IR limit of the broken phase independently of the RG scheme. In the symmetric phase the 
evolution of z is negligible in the IR giving the same evolution as was obtained in LPA with the line of fixed points. 

One can easily show that the critical exponent r/ v characterizing the vortex-vortex correlation function [144] is 
rj v = 1/4 independently of the parameter b [28] . However the anomalous dimension for the correlation function of 
the field variables gives i] = in the vicinity of the KT point. In the deep IR scaling region the situation changes 
significantly, there new scaling laws appear. Fig. 14 shows that around the KT point (at about fc/A ~ 1CP 4 ) z is 
practically constant, giving r\ — 0, while in the IR region z starts to diverge according to a power law scaling. This 
scaling regime is induced by the IR fixed point. The corresponding exponent depends on the scheme, i.e. on the value 
of b. The coupling u shows the same qualitative behavior. If one plots z as the function of k then one obtains the 
results of Fig. 15. The flows show sharp singularities as the function of the scale k so as the other couplings, making 
the evolution scheme independent. The modes which are integrated out during the RG procedure can be indexed not 
only by the scale k but the scale k and one has a smallest value of the scale k c . Again one can identify the correlation 
length according to £ = l/k c , furthermore the critical exponent also can be identified in a similar manner as was 
introduced in the d-dimensional O(N) model. In the 2d SG model the temperature is proportional to Zfe A , and thus 
the reduced temperature t is t ~ Zk A — z c k A , where z c & A is the initial value of the wavefunction renormalization at 
the separatrix. In order to demonstrate whether the phase structure is qualitatively the same we plotted it for b = 5 
in Fig. 16. In its inset we plotted log£ as the function of t. We obtained a straight line implying the scaling of £ 
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k/k A 

FIG. 14: The scaling of the wavefunction renormalization z is plotted as the function of k. The initial values are: z{ktC) = 0.6 
and u{k\) around 1.95. 




k/k A 

FIG. 15: The scaling of the wavefunction renormalization z is plotted as the function of k. The initial values are: z{k\) — 0.6 
and u{kh) around 1.95. 



according to 

log^cx* - ". (76) 

There are two types of correlation lengths, one is defined as usual, i.e. around the KT turning point of the coupling 
u. Another one is identified as £ = l/k c in the neighborhood of the IR fixed point. It can be seen from the inset of 
Fig. 16 that the scaling of £ shows an infinite order phase transition for all schemes with the exponent v w 0.5. 

We obtained that the exponent v coincides calculated around the KT and IR points. This is due to fact that 
condensate is a global object in the broken phase which is perceptible at practically any scale. The other exponents 
do not necessarily coincide at the KT and the IR points, e.g. the anomalous dimension -q = around the KT point, 
while it is t] = 2b/ (b — 1) around the IR one. 

The singularity at low scales appears because the condition in Eq. (75) is satisfied. The condition shows that 
there are many modes in the model with infitesimally small, practically zero energies, they are called soft modes. 
The ground state of the model can pick up arbitrary number of these fluctuations, which makes the ground state 
degenerate. The huge amount of soft modes constitute the condensate in the broken phase. They signal that the 
microscopic degrees of freedom are not suitable to describe the model anymore in the broken phase at those small 
energy scales, one should turn to new field variables. It is reminiscent of the composition of hadrons from quarks, 
which are the original degrees of freedom in QCD. 

One can conclude that the IR fixed point provides us the low energy limit of the theory under investigation. On 
the other hand the soft modes suggests that the quantum-classical transition appears in the broken phase, however 
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FIG. 16: Phase diagram of the SG model, with b — 5. The dashed (solid) lines represent the trajectories belonging to the 
(broken) symmetric phase, respectively. The wide line denotes the separatrix between the phases. The KT point can be 
considered as the crossover fixed point. The inset shows the scaling of the correlation length as the function of the reduced 
temperature t. The curves are shifted for better visibility. The lower (upper) set of lines corresponds to the IR (KT) fixed 
point. The triangle, circle and square correspond to b = 2, 5, 10, respectively. In the middle a straight line with the slope —1/2 
is drawn to guide the eye. 



correct treatment should be discussed by the close time path method [17, 145]. 

UV NGFP 

The UV NGFP also cannot be seen directly analytically but it can be easily recognized in the top right corner of 
the phase diagram in Fig. 13. The LPA approximation trivially shows that in the UV attractive NGFP the model 
has a singularity, as in the IR fixed point. The UV singularity signals the upper limit of the applicability of the 2d 
SG model. It implies that at high energies the model requires new elementary excitations. One can make it more 
apparent in the framework of the XY model, which is in the same universality class as the 2d SG model. There the 
excitations are represented by vortices consisting of concentric forms of spins [146]. The blocking towards the higher 
scales means that the vortices have smaller and smaller diameter. Naturally at a certain scale we reach a limit, where 
the vortex reduces to a single spin. At this scale the new elementary excitations should be the single spins instead of 
the original vortices. Therefore in the UV limit the system of the charged vortices should be replaced by a neutral 
spin system. 

The 2d SG model shows a nice example, where the functional RG method gives us both the high and low energy 
scale limits of its applicability. The low energy limit is usually indicated by the IR fixed point as was demonstrated in 
the article for many models. There the global interaction belonging to the condensate of the broken phase becomes 
IR relevant and introduces new excitations. However there should be an UV limit of the models, since it seems 
nonphysical to assume that a model can be valid at arbitrarily high energies. So far the functional RG method gave 
the upper UV limit only for the 2d SG model and its generalizations (massive [147, 148] or layered [86, 149-153] SG 
models). 

D. Quantum Einstein gravity 

The QEG is the quantum field theoretic model of gravity [93, 107]. There the metric fields are considered as 
fundamental degrees of freedom, and play the role of the field variable, thus the model can be treated in the framework 
of path integral formalism. The simplest version of the model contains the Newton constant Gk and the cosmological 
constant A& as couplings. According to their canonical dimensions the dimensionless cosmological constant A = k~ 2 Ak 
is relevant, i.e. perturbatively renormalizable, however the dimensionless Newton constant g scales in irrelevant 
manner, since g — Gk d ~ 2 . Thus the model seems to be non-renormalizable according to perturbative considerations. 
This might be the signal of the necessity to introduce new dynamical variables for elementary excitations, e.g. strings. 

However the concept of asymptotic safety may give a helping hand to QEG. It was shown that the model has an UV 
attractive NGFP, which makes QEG renormalizable [78, 88, 89, 91, 96] and shows asymptotic safety [103, 104, 154- 
158]. It is a great challenge to find experimental evidence for quantum gravity, however there are promising ideas to 
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catch its effects in LHC [159, 160]. Usually the model is considered in d — 4, however the 3-dimensional [161] and 
higher dimensional cases [162-164] are also thoroughly investigated. In the context of dimension we notice that a 
dimensional reduction takes place during the flow to small distances from 4 — > 2 [93, 165-167]. The bare action of 
QEG belongs to the UV attractive NGFP of the RG flows [93, 167]. Some kind of classical limit can be associated to 
any fixed points, which can give different General Relativity and cosmology [157]. 

QEG is defined through a diffeomorphism invariant functional of the metric g^ v in QEG, and the background 
field technique is used to preserve the gauge symmetry when the RG flow equations are derived. Here we flash the 
main points of the derivation given in [78]. The functional integral should be performed over all metrics 7 M „. It is 
decomposed into 

l^v(x) = g^ v (x) + h^ v (x) (77) 
where g p „{x) is the background field metric and h^ u {x) is the new integration variable. The generating functional is 

MW k } = J V[h^;C;C}e X p{S\S + h]-S gf [h;g}- S gh [h,C,C;3] ~ &kS[h,C,C;g] - S source }, (78) 



exr 



where S[g + h] = S[j] is the classical action, which is invariant under the general coordinate transformation (C v is 
the Lie derivative w.r.t. the vector field v) 

<hv = £v = v p d p jf, u + d p v p ^ pu + d u v p jn P . (79) 

The term S g f in Eq. (78) denotes the gauge fixing term 

S gf [h;g] = ^J d d xV§g^F^F u , (80) 

where F p is linear in the field variable h^ v and contains a certain form of a first order differential operator of g pu [78] . 
The term S g h stands for the Faddeev-Popov term with the ghosts C M and 

S gh [h,C,C;g] = — / ^xC^j^Ccig^ + h^), (81) 

with n — (327rG)~ 1//2 . It is obtained by applying the gauge transformation 

Shpv = CvJuv, Sg^ = (82) 

to F p and then replacing the parameters w M by the ghost field C*. The Faddeev-Popov determinant is represented 
by the path integrals over C M and C p . The IR regulator should be applied to both the ghost and the gravitational 
field, and it reads as 

A k S[h,C,C;g] = ±k 2 J d d x^§h^K° rav [gr P °h pa + V2 J d d x^C^ h [g]C^. (83) 
The source terms in Eq. (78) has the form 

d d xV3(t^h^ + + o»Cp + P^Ccifa + V) + r p C v d v C»). (84) 



In order to get the effective action one should introduce the classical field variables 

- J_5Wk_ _ J_SWk - J_SWk 
M " ~ ^g 6t»» ' ^ ~ sfg da"' ^ ~ VI ^ ' 

and then it can be obtained by the following Legendre transformation 

r fe = j d d xV§(t^h^ + + tr%) ~W k + AS k . (86) 
Finally one arrives at the flow equation of the form 



t*\h £ £ Si = -Tr ^ n9 k aV ^h i T 



> 2 ^f +no« (* 2 n h +n\L 



(87) 
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A general form of the QEG effective action is 

T k = [ d*Xy/d^Z (-nT7^( 2A k - Ft) - + ■■■ + t^~C 2 + —£+... + gf. terms + gh. terms) . (88) 

J \167rG fe 6a k la k a k J 

The first terms come from the Taylor expansion of the general non-local functional f(R) [168, 169]. The first two terms 
constitute the Einstein-Hilbert truncation, with the two dimensionless couplings A = A k k~ 2 and g — G k k d ~ 2 . Then 
the contributions coming from the square of the Weyl tensor C 2 = C ilvpa C >lvpa ', and the Gauss-Bonnet topological 
invariant E = R^ pa -R fJ ' uf " T — 4i? M ^i? AlI/ + R 2 are considered [170] in Eq. (88). In the end the gauge fixing and the ghost 
terms arc indicated. The functional f(R) can have any form, e.g. it can be a logarithmic or polynomial function 
[168]. The explicit form of the evolution equations including the interaction term R 2 can be found in [171]. A further 
possible extension of QEG contains terms in the action in Eq. (88) which describes interaction between the matter 
and the metric fields. The matter field can be a fermionic system [172-174] a scalar field [102, 175-180], or a gauge 
field [108, 181, 182]. The improvement of the ghost sector is also intensively studied [183, 184]. 

From the extensions of the QEG effective action in Eq. (88) one can derive the RG evolution equations for the 
couplings. Considering the 2-dimensional projection of the phase space to the usual Newton and cosmological couplings 
we can obtain different types of phase structures. However it has several universal properties, namely that any 
extension of QEG shows two phases and a possesses a UV attractive NGFP, implying that model is asymptotically 
safe. 

E. Evolution equation 

Using the Litim's regulator and the Einstein-Hilbert truncation, where the dimensionless Newton and the cosmo- 
logical couplings are considered, we obtain the following RG equations [164] 



A = -2X+"-d{d + 2){d + b)-d(d + 2Y- 



l-4A(l-l/d) 
d-2 



2 y " ' v '2 g-g b 
d + 2)g 
9- 9b 



with 



One can introduce the gravitation anomalous dimension 

„ = (91) 

9 - 9b 

The model has two fixed points [162, 163, 185]. When d = 4, then there is a UV NGFP at \* uv = 1/4, g* uv = 1/64. 
The matrix M is 

«4(r-7)' 

with the eigenvalues Suvi — ( — 5 + i\/167)/3 and sjjv2 = (—5 — i\/l67)/3, so it is an UV attractive of IR repulsive 
focal point. Furthermore one has a GFP at Xq — 0, g\, v — with the corresponding matrix 

M= l T h W + W-*)\ (93) 
3 \ — 1 d — 2 J 

in dimension d. The eigenvalues are sqi = —2 and sq2 = d — 2. The negative reciprocal of the former eigenvalue 
gives the critical exponent v = 1/2 of £. The latter one guarantees that the GFP is a hyperbolic one when d > 2. 

The phase structure is shown in Fig. 17. There is a separatrix connecting the UV attractive NGFP and the GFP. 
The trajectories which are left to the separatrix give negative values for the cosmological constant and vanishing 
Newton coupling in the IR limit, and constitute the strong-coupling or symmetric phase of QEG [78, 186]. Other 
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X 

FIG. 17: The phase structure of quantum Einstein gravity in d — 4 is shown. There is an UV attractive NGFP, a crossover 
fixed point, which is now the Gaussian one, and an IR fixed point. The thick line represents the separatrix. 




K -K 

FIG. 18: The scaling of the correlation length, giving v = 1/2 for the exponent. For simplicity we fixed the UV value of g\. 



trajectories getting around the separatrix from the right give large positive values of A and small Newton coupling 
when the RG flows approach the IR regime. This phase is called the weak-coupling or broken symmetric phase. 

In d = 4 if one reparametrizes the couplings according to x = 1 — 2A, to = 4g — (1 — 2A) 2 and introduces the new 
'RG time' d T — ud t then one obtains 

dr X - -4^ + 2x^(8 + 21 X ) + 24 W 2 + 6x 2 (3 X (x+l)-l), 

d T u = 8^ 2 (1 - 6 X ) - 2 X (42 X 2 + 9 X - 4) - 6 X 3 (x(6x + 5) - 2), (94) 

which have the UV NGFP at oj'(j V = —3/16 and xhv = V^> the GFP at ojq = —1 and xh = 1 an< ^ the IR attractive 
fixed point at ui} R = and x*ir = which corresponds to the point g* IR = and A}^ = 1/2 in terms of the original 
dimensionless couplings. Naturally the IR fixed point is IR attractive. The existence of the IR fixed point also has 
been uncovered and discussed in [108-110]. 

As in the case of the previous models, one can get the value of the critical exponent v from the IR scaling. Again, 
the scale k where the evolution stops equals the reciprocal of the correlation length £. We plotted the divergence 
of the correlation length in Fig. 18 for dimensions d = 4... 11. We obtain that a second order phase transition 
can be identified in the IR region of QED with the exponent v = 1/2 independently on the dimension. One can 
obtain similar results if one uses different regularization scheme or special extensions of the action containing the 
functions of the Euclidean spacetime volume V = J d d x Actg ^ v [87]. This result does not coincide to the UV value 
of v = 1/3 [187, 188], which is the reciprocal of the imaginary part of the scaling exponent in the UV fixed point, 
i.e. v = 3/\/167 « 0.23. They are not equal but taking into account further couplings e.g. by including higher order 
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model 


UV 


CO 


IR 


3d O(N) model 


Gaussian 


Wilson-Fisher 


IR 


3d nonlinear a model 


non-Gaussian 


non-Gaussian 


Gaussian 


3d Gross-Neveu model 


Gaussian 


non- Gaussian 


IR 


2d sine- Gordon model 


Gaussian and non-Gaussian 


Kosterlitz-Thouless 


Gaussian and IR 


4d quantum Einstein gravity 


non- Gaussian 


Gaussian 


IR 



TABLE III: Summary of the models and their fixed points. 



terms in the curvature scalar [118, 189-194] the result of v can be improved. Generally the exponent has some scheme 
and truncation dependence, and it is around v w 0.23 — 0.4. 

However the improvement of the action by extensions does not necessarily change the IR value of v = 1/2. As it was 
demonstrated through many models in the previous sections the IR value of v equals the one which can be calculated 
in the vicinity of the crossover hyperbolic fixed point if it exists. Here the GFP plays this role, but we know that 
the eigenvalues of the GFP equals the negative canonical dimension of the couplings. The single negative eigenvalue 
of M in Eq. (93) is s\ — —2 independently on the dimension, as is numerically checked and shown in Fig. 18. If we 
introduce new couplings into the QEG effective action in Eq. (88) as a Taylor expansion of the dimensionless form of 
f(R) [195] according to 

W = ES^ ( 95 ) 

n=l 

then the canonical dimension of the coupling g n becomes d n = d — 2n. We can identify the first two dimensionless 
couplings with the cosmological and the Newton couplings as A ~ go jg\ and g ~ 1/fifi • Then the canonical dimension 
induced scaling of the couplings around the GFP are A ~ fc -2 , g <~ k d ~ 2 and g n ~ k 2n ~ d with n > 1. In case of d = 4 
the single negative eigenvalue is s\ = —2, which belongs to the cosmological coupling. It implies that the exponent v 
calculated at the GFP is always equals 1/2, therefore we expect, that its IR value is exactly v — 1/2 if the total f(R) 
is considered. 

We note that there are extensions of QEG where the GFP does not exist [107]. Since the IR fixed point is strongly 
related to the crossover hyperbolic point then its scaling behavior can significantly change. It was shown that a 
certain choice of the extension may eliminate the GFP [87]. In this case the correlation length £ does not diverge 
which implies that the order of the phase transition changes to a first order one, in the IR region. Naturally it can 
also change the classical limit there which may also affect the cosmological value problem [87]. 

IV. SUMMARY 

We gave a short introduction of asymptotic safety in the framework of the functional renormalization group method. 
We introduced several models, where a nontrivial UV attractive fixed point exists, which saves the model from 
nonphysical divergences and gives asymptotic safety. 

The functional renormalization group method is a powerful nonperturbative tool to investigate the field theoretical 
models. It provides a partial differential equation for the effective action which initial condition describes the high 
energy UV physics of the model and the solution gives the low energy IR one. The RG technique is used in several 
areas in quantum physics, however the results from its nonperturbative nature comes forward quite rarely. One can 
say that the stable UV and the stable IR limit of the theory cannot be described by perturbative calculations. The 
former one provides the asymptotic safety, while the latter one gives the IR behavior of the broken phase which can 
be related to the quantum-classical transition. In this articles we concentrated on these limits of the models which 
were considered. 

First we investigated the d-dimensional O(N) model, which is asymptotically free. We showed that there exists two 
phases and three fixed points in the model. Then we analyzed such models which are asymptotically safe. We plotted 
the phase structure of the 3d nonlinear a, the 3d Gross-Neveu, the 2d sine-Gordon models and the 4d quantum Einstein 
gravity. These models have two phases. The phase space shows several similarities. There are three different fixed 
points. The trajectories start from the UV attractive NGFP and tend towards the crossover hyperbolic fixed point. 
In the broken phase of the models there exists a further IR fixed point which is IR attractive. Table III summarizes 
the fixed points of the models, which were investigated. The models are not in the same universality class, since they 
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can be distinguished by the critical exponents. Different fixed points may result in different exponents, therefore the 
classification could be done for the UV, CO and the IR fixed points. 
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